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Abstract 



A consistent perturbation theory expansion is presented for phase ordering 
kinetics in the case of a nonconserved scalar order parameter. At lowest order 
in this formal expansion one obtains the theory due to Ohta, Jasnow and 
Kawasaki (OJK). At next order, worked out explicitly in d dimensions, one 
has small corrections to the OJK result for the nonequilibrium exponent A 
and the introduction of a new exponent v governing the algebraic component 
of the decay of the order parameter scaling function at large scaled distances. 
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Typeset using REVT^ 



Significant progress has been made on the theory of phase ordering kineticsS using meth- 
ods that introduce auxiliary fields that are taken to have gaussian statistics. The methods 
developed in OJkI and TUci each have appealing aspects and separately give good de- 
scriptions of different aspects of the ordering problem. A major lingering question is why 
these methods work as well as they do and how they can be reconciled and improved. Thus 
there has been a searchSi for the field theory description where OJK or TUG is the zeroth 
order approximation in some systematic expansion. Such an expansion is presented in this 
paper. It has the OJK result as its zeroth order approximation. 

The system studied here is the domain wall dynamics generated by the time-dependent 
Ginzburg-Landau model satisfied by a nonconserved scalar order parameter ilj{f,t): 

A(i)z^(i) = -v'mi)] (1) 

where, in dimensionless units, A(l) = ^ — is the diffusion operator and 1 denotes (ri, ti). 
We expect our main results to hold for any symmetric degenerate double- well potential V[ip] 
and to be independent of the exact nature of the initial and final states, provided they are 
disordered and ordered states respectively. Therefore it is convenient to quench to zero 
temperature and set any noise term to zero. 

It is well established^ that for late times following a quench from the disordered to 
the ordered phase the order-parameter correlation function has a scaling form C(12) = 
(■?/'(l)-?/'(2)) = iPqJ^{x, ti/t2) where ipo is the magnitude of the order-parameter in the ordered 



phase and the scaled length x is defined as x = {fi — f2)/L where L = v4T = y 2(ti + 
is a growing characteristic length. In the perturbation theory developed here we find that 

^(x, h/t2) = - sin-^ [ $(ti/t2)/(a;, h/t^)] (2) 
vr 

with the normalizations /(0,ti/t2) = 1 and $(1) = 1. The quantity $(^1/^2), related to the 
on-site order parameter autocorrelation function, has the formi 

^itilt2)=i^\ (l + A$(ti/t2)) (3) 



where A$(ti/t2) is a regular function of ti/t2 which is of second orderi in the expansion. 
The nonequihbriumi exponent A, which gives the decay for ti ^ t2 and vice versa, is given 
to second order by 

The dimensionahty dependent quantities, tu, and M^, are given by 

^ f Md \ d 



and 



2u; + u^2^K,+ ^^)=l + - (5) 



Mrf = / dz- -T = Vir^ • (7) 

fix^txjt'}) can be written in the form 

where A/(x, tijt-i) is a regular function of its arguments and is second order in the expansion 
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and /o(a;) = e~ 2^' is just the O JK result for the scaling function. The exponent v governing 
the algebraic component! of the large x behavior, shared by / and JF, is given by 

. = + . (9) 

At lowest order one has the O JK results A = | and z/ = 0. While Kfi and uj can be worked out 
analytically for specific values of rf, the expressions are not very illuminating. Numerically 
we have for = 1, = 0.4601.., A = 0.6268.., v = 1.1596..; d = 2, u = 0.6877.., A = 1.1051.., 
u = 1.2492..; d = 3, cu = 0.9067.., A = 1.5828.., u = 1.3732... For large d, to ^ 1(^2-1), 
A ^ f, and z/ ^ d{3-2V2). The 6est estimates^ are ci = 1, A = 1, = 2, A = 1.246 ±0.02, 
and d = 3, X = 1.838 ± 0.2. For equal times ti = t2, T{x^ 1) = T{x) = | sin~^f{x, 1). 
For small x, where / is analytic in an expansion in even powers of x, satisfies Porod's 
law0 and the Tomita@ sum rule. 



The perturbative treatment developed here is a two step process. The first step is to 
introduce an auxihary field m(l) which satisfies an equation of motion 

A(l)V,m(l) = 2i{h)pi{l) = 2i{tl)5{m{l))V^m{l) . (10) 

The quantity ^(ti) must be chosen to go as L^^{ti) for large times with a coefficient which 
must be adjusted if the system is to order, much like the temperature must be set to in 
critical phenomena. This equation of motion has the appealing feature that Vim(l) is a 
diffusion field driven by a source pj(l) which is the invariant density^ of domain walls for 
the system. We show below that we can organize a perturbation theory expansion for this 
field theory about a gaussian zeroth order theory. 

The second step in the analysis is to show that one can map the phase ordering kinetics 
problem onto this new field theory. This mapping is via the transformation if) = a[m] + u[m] 
where (j[m] satisfies the Euler-Lagrange equation^ for the associated stationary interface 
problem with m the coordinate: cr2 = V^'[c'"[m]] where we introduce the notation ai = 
The idea is that m grows large away from an interface and a —>■ tpQsgn{m) in the bulk. 
Inserting this mapping into the equation of motion for ip we find: 

A(l)M(l) + ai(l)A(l)m(l) = -l^'[a(l) + n(l)] + a2(l)(Vm(l))2 . (11) 

The most important property of the solution of this equation for u[m] is that Zzm|m|-^oow[m] = 
and — > ipQsgn m. In the regime where |m| is large, o"^ ~ tp^, the derivatives of cr go 
exponentially to zero, the linearized equation for u shows that it has acquired a mass and, in 
the long-time long-distance limit, u is seen to go exponentially to zero. The universal prop- 
erties of interest are associated with bulk quantities like C{12...n) =< iIj{1)iIj{2) ■ ■ ■ip{n) > 
where the points ri, r2, r„ are not constrained to be close together. Because the field u{l) 
is nonzero only near interfaces the average < u{l)a{2) ■ ■ ■ > is down by a factor of 1/L^(ti) 
relative to < a{l)a{2) ■ ■ ■ >. Thus if we focus on these properties we do not need to know 
u in detail. 

We will develop the field theory generated by the equation of motion, Eq.(|r^), for m(l) 
using the standard Martin-Siggia-RoseEl method in its functional integral form as formulated 
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by DeDominicis and Pelitilli. In the MSR method the field theoretical development requires 
a doubling of fields to include a field rhi conjugate to V^m. We also organize things so 
that the initial field mo(r) is also treated as independent. Following standard procedures, 
averages of interest are given as functional integrals over the fields m, rhi and mo weighted 
by an action A: 

A{m, m, mo)=i j dlY, m,(l) [A(l) V,m(l) - 2^(1)^(1)] (12) 

i=i 

-i J c/l^mi(l)(5(ti -to)Vimo(l) - ^J d\i J mo{fi)g-\fi - r2)mo(f2) 

i=l 

where we use the notation dl = dtid'^ri and where are assume, as is appropriate in this case, 
that the initial field is gaussian and has a variance given by < m,o(r*i)mo(r*2) >= di'f^i — ^2) 
. We will not have to be very specific about the form of the initial correlation function g. 
It is convenient to generate our correlation functions as functional derivatives in terms of 
sources, h{l) and h{l), which couple to the conjugate fields. Thus we introduce the total 
action At'- 



At = A + j dl 



h{l)m{l) + h{l)Y,Vim,{l) 



i=l 



(13) 



The fundamental equations of motion are given then by 



< '^^J^^Tirn, m, mo) >h= (14) 

where ^1/(1) = {m(l), mj(l), mo(l)} and the averages are over the total action At- Taking 
the functional derivatives of At and and, using Eq.(0) with \1/(1) = mo(l) to eliminate the 
initial field in terms of the hatted field, we obtain the two basic equations 

A(l) < M(l) >h +2^(1) < p{l) >h] = h{l) (15) 

I [A(l)Vi < m(l) >h -2^(1) < Pi{l) >h] = -Vf ^ I d2 no(12) < M(2) >h +vMl)- (16) 

Here we have introduced the quantities no(12) = 6{ti—to)5{ti—t2)g{ri—f2) ,A(1) = ^+Vf , 
p{l) = 5(m(l))M(l), and M(l) = Vimj(l). Notice that the fundamental equations 



Eqs.([T5|) and (|^) depend on M and not rhi and we have constructed h such that it couples 
to M. All correlation functions of interest can be generated as functional derivatives of 
< m(l) >h or < M(l) >h with respect to h and h. 

Taking the functional derivative of Eqs.(|l^) and (0) with respect to h{2) gives the 
equations for the two point response and correlation functions: 



A(1)G^^(12)+Q2(12) =5(12) 



(17) 



A(1)vS')G^„(12) -Q^^^(12) = -vr / dS no(13)G^,„(32) 



(1) 



where we have used 

Qn(12...n) 
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n-l 



6h{n)6h{n-l)---6h{2) 



[2^(1) <p(l) 



Q[:ni2. 
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n-l 



.n] 



[2^(1) < a(1) >,] 



6h{n)6h{n-l)---6h{2) 
evaluated for n = 2 and introduced the notation 

< M(l) >^ 



<^Mmm..m(123..n) 



5hi2)5h{3) ■ ■ ■ 5h{n) 



GJl2---n) 
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n-l 



< mil) >, 



(18) 



6h{n)6h{n - 1) ■ ■ ■ 6h{2) 
The functional derivative of Eq . ([T6|) with respect to h{2) is redundant compared with Eq.(|lj 
because of the relation G'^jyj(12) = Gj^j,^{21) . The equations governing the n^'^ order 
cumulants are given by 



A(1)G 



Mm...m 



;i2...n) +g„(12...n) 







A(l)V.G„(12...n) - Q»(12...n)] = -VS^^ / dl no(ll)G^,^...^(12...n) . 

The point now is to show that there is a consistent perturbation expansion for this theory. 
After introducing the one-point probability distribution Ph{xi, 1) =< 6{xi — m(l)) >h, it is 
easy to show 



6 



Q?{1) = 2^(1) < a(1) >h= I dxia^)sgn{xi)Ph{x,, 1) 



(1) 



(19) 



gi(l)= 2^(1) </)(!) >h= 2^(1) 



< M(l) >h + 



PhiO,l] 



(20) 



5h{l)\ 

Then any perturbation theory expansion for Ph{xi, 1) will lead immediately to an expansion 
for Q^^ and Qn by functional differentiation. The expansion for P/i(xi, 1) is straightforward. 
Using the integral representation for the 5-function we have 



>h-- 



dki 



e ^ ^exp 



.n=l 



n 



27r J 27r 

where H{1) = ikim{l) and G^J}\l) is the n*^ order cumulant for the field H{1). Since H{1) 
is proportional to m(l) these are, up to factors of iki to the n^^ power, just the cumulants 
for the m field. It turns out that in zero external field we can take Gn to be of (9(|^ — 1). 
Keeping terms up to the 4-point cumulant it is easy to see that 



(21) 



where 



Pi'\x„l) = J^Mh,l)e 



-ik\x\ 



(22) 



and 



$o(fci,l) = e 



-ifci(xi-Gi(l))g-ifc2G2(ii) 



(23) 



Eqs.([2l|), (p^) and (^) inserted back into Eqs.(|l9]) and (pOD gives a relation between the 
m-field cumulants and the quantities Q^^ and Qn- These quantities are ultimately related 
to derivatives of -P/[°''[0, 1]. Because of this we see that each term in the perturbation theory 
expansion in zero external fields is weighted by factors of 



0n(l) 



dki 
27 



d 



'11 



(24) 



where 



1) = ^ s{i) ^"^^ ^^^^ m^(l) >. The structure of the perturbation theory 
becomes clear if one looks at two classes of contributions to Qgn- The first class corresponds 



to the 2n — 1 derivatives which all act on $o(^i; 1) to give contributions to Q2n(12 ■ ■ ■ , 2n) 
of the form 

vS')e(l)0n-i(l)G'2(12)G'2(13) ■ • ■ G2(l, 2n) . 

There is another set of contributions where all of the derivatives except the first acts on the 
factor multiplying $o(^i) 1) and gives a contribution to (52n(12 ■ ■ ■ , 2n) of the form 

vS')e(l)0o(l)G'2„(12---,2n) . 

Inserting these results into the expression for the 2n point cumulant we see term by term 
that the leading order behavior goes as ~ 4'n-i- Orders in the expansion correspond to 
the sum of the labels n on 0„ in products of 0's. Thus a term with factors 0i020i, typically 
associated with different times, is of 0{A). It should be emphasized that at this stage that 
this is a /orma/ expansion. At order n it is true that (pn ~ L~'^^"+^) is small, however it will 
be multiplied , depending on the quantity expanded, by positive factors of L(t) such that 
each term in the expansion in 0„ has the same overall leading power with respect to L{t). 

It is then straightforward^ to work out the expansion for the two-point cumulant and 
obtain at zeroth order: 

g2(12)=^o(l)G'^,™(12) (25) 

Q«(12)=VS%(1)G2(12) (26) 

where we have defined co'n(l) = 2^(l)0n(l)- The equations of motion at zeroth order for the 
two point response and correlation functions is given by Eqs. (|1^ and (|18D with Q2 and 
given by Eqs. (|25|) and (pBI): 

A{1) + Ml)] G^^JU) =6(12) . (27) 



and 



t [A(l) - ^o(l)] (12) = - / dl no(ll)G2^^(12) . (28) 



The equation for the response function has the solution 



where we have introduced 

Rihh) = eK '^^''^^^ . (30) 

The equation for the correlation function, Eq.(^), can be integrated up into the symmetric 
form, suppressing the step functions in time, 

Gf\r,hh) = Rihto)R{hto) J ^,e^'-%q)e-^''"^ (31) 

where T = and g{q) is the Fourier transform of g{r). We focus on the long-time limit 
where in the integral over wavenumber we can replace g{q) g{0) to leading order in powers 
of 1/T and obtain for large T: 

Gf\r,hh) = R{tito)R{t2to)m ■ (32) 

If this system is to order, < m^(l) > is to grow large, then we must choose, for large time 
t, uJo{t) = 2^(t)0o(t) = J where a; is a constantS we will determine. Then, R{tit2) defined 
by Eq.(pO|), is given by R{Ut2) = ('f-T and 



If we are to have a self-consistent scaling equation then the autocorrelation function (r = 0), 
at large equal times ti = t2 = t must satisfy 

Gf (0.«) = S<»>(i) = A„t = i--/^^^ (34) 

which fixes the exponent oj = |(1 + |) and the amplitude Aq. This determination of 
uj must be repeated order by order in perturbation theory. Returning to the correlation 
function, we can eliminate ^(0) using the expression in terms of the amplitude and obtain 

G^2\r,hh) = ^S{h)Sit2)%ih/t2)fo{x) where $0(^1/^2) = (4^)'^' and /o(x) = e^^'/^. 
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One can read off from these results the nonequihbrium exponent A = | given by the OJK 
result. 

Going forward, it is easy to find that the next nonzero contributions to Q2 and Q2 is of 
0(2) and given by 

(1) r^i(i), 



3! 



-^4(1112) 



(35) 



4(12) 



(2) 



^l(l) 



^Mmmm(1112) 



(36) 



These require that we evaluate the four point cumulants at the leading first order and obtain 



GM_(1234) = G£)^(li)Ha;i(i 



(37) 



X 



G'2jl2)Gf\l2>)Gf\U) + Gll(13)Gf (i2)Gf (14) + G^^l{U)Gf\l2)Gf\U) 



Mm 



6-4(1234) = G'f(li)Hu;i(l)] 



G^SJ12)gP{13)G^^\u) 



(38) 



+ GS^^(i3)Gf (i2)Gf (14) + G^°^^(i4)Gf (i2)Gf (13) 



+ ^!!!w(ll)(-^^i(l))^2°Hi2)G(°)(i3)G?)(i4) 

where we integrate over the repeated barred index. Notice that G4 is properly symmetric 
under interchange of its labels. Inserting these results back into the equations for the two 
point quantities we have for the response function 



GmJ^2) = G^^JU) + / dld2G£)^(ll)E^^(12)G£)^(22) 



amd the lowest order self-energy contribution is given by 



^Mmm = §Hc.i(l)]GSL(12)Gr(12)1-ia;i(2)] 



.(0)/ 



The correlation function is then given to second order by 



(39) 



(40) 
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G2(12) = (12) + G^''')(12) + ^^'''^(21) + ^^'''^(12) (41) 



where 



G?'^)(12) = / rfld2Gf (11)S,^„(12)G(°)^(22) (42) 



and 



G?^)(12) = -/ rfid2G;;). (li)n(^)(12)Gg)^(22) (43) 

where the self-energy is the same as for the response function and 

n(2)(12) = -lH^i(l)]G'f(12)=^H^i(2)] . (44) 

The spatial integrals giving Gf^^ and G^^'^^ can be evaluated in d-dimensions since they 
involve products of displaced gaussians. After rescaling the internal time integrations ii = 
Tyi and ^2 = Ty2 we obtain 



Gf'')(12) = ^5(1)^(2)2^- V$o(ti As) Ji(a;, ti/T, T) (45) 

ito/T Jto/T [(yi + i/2)(3i/i - 1/2 - (2/1 - y2r)r 



d-l 



Gf''\l2) = ^S{l)S{2)'^uj^%{ti/t2)J2{x, h,t2) (46) 

fti/T rt2/T y'l^'^^^y^'^^^ 1 3 ^2 

Ao/T Ao/T [(?/i + y-iYKi - 2/1 - 2/2)]^^/^ 

The first thing we should do with these results is look at the contribution at this order to 

the on site equal-time ti = t2 = t correlation function. The integrals Ji(0, 1, t) and J2(0, t, t) 

are logarithmically divergent as t ^ 00 and we have 

^(1) = ^W(t) [1 + (k, + ln{tlt,) + ■ ■ -l (47) 



3d/2 

where Kd is defined by Eq.(|^) and Md is defined by Eq.(|^. Remembering that at lowest 
order S^^\t) = AqI^^''^/'^ and after exponentiation of the ln{t/to) term in Eq.(^7D we find 
that to obtain S{t) = At, we require that uj be given by Eq.(|^). 
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The to singularities in Ji and J2 can be regulated by turning our attention from 
6*2(12) to the scaled quantities $(^1/^2) and f{x,ti/t2) which appear in Eq.(^: 

/(12) = ^S2|= = $(ti/t2)/(x,ti/t2) . (48) 

Looking at the second order contributions to $(^1/^2) we find for ti ^ t2 or t2 ^ ti there 
are contributions which are logarithmically divergent in the quantity ^/t^jT which can 
be exponentiated to obtain the result given by Eq.(^. f{x,ti/t2) is an analytic function 
of X. One can easily work out, for example, the lowest order term in the power series 
expansion in for f{x,ti/t2) which is regular for all d. It is this term which gives the 
coefficient in Porod's lawli. There is one last regularization that must be carried out before 
the perturbation theory expression for the correlation function can be used for all values of 
X. If we look at large x we find that f{x,ti/t2) is logarithmically divergent in 1 + which, 
on exponentiation leads to the expression given by Eq.(^. 

We turn next to the connection between the correlation function for the auxiliary field 
m and the order parameter correlation function. As discussed earlier, since u{l) vanishes 
exponentially for large |m(l)|, then the averages over these fields are down by a factor of 
relative to the averages over the field a and in the scaling regime C(12) = {a{l)a{2)). The 
perturbation expansion for this quantity follows closely the development for the one-point 
quantity except we must treat the two- point generalization of Ph{xi, 1). One can then show 

C(12) = - sin-^ /(12) + 0(3) (49) 

where /(12) is defined by Eq.(^). Since we have worked out /(12) to second order previ- 
ously, we have consistently to this order, the result given by Eq.(|D. 

Clearly one can go to higher order in this perturbation theory and introduce more sophis- 
ticated graphical methods. It is interesting to ponder the implications of the ideas developed 
here for the conserved and vector order parameter cases. 
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